Mixed-phase description of colossal magnetoresistive manganites 
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In view of recent experiments, indicating the spatial coexistence of conducting and insulating 
regions in the ferromagnetic metallic phase of doped manganites, we propose a refined mixed- 
phase description. The model is based on the competition of a double-exchange driven metallic 
component and a polaronic insulating component, whose volume fractions and carrier concentrations 
are determined self-consistently by requiring equal pressure and chemical potential. The resulting 
phase diagram as well as the order of the phase transition are in very good agreement with measured 
data. In addition, modelling the resistivity of the mixed, percolative phase by a random resistor 
network, we obtain a pronounced negative magnetoresistance in the vicinity of the Curie temperature 
To. 

PACS numbers: 71.10.-w, 75.47.Gk, 71.38.Ht, 71.30.+h 



I. MOTIVATION 

The peculiar properties of the ferromagnetic metallic 
phase of mixed- valence manganites (e.g. Lai-^Ca^MnOs 
with 0.2 < x < 0.5), in particular the large negative mag- 
netoresistance close to the Curie temperature Tc, have 
been the subject of intense research activity over the last 
years 1,2,3,4 . There is now compelling experimental ev- 
idence, that the complex interplay of the electronic de- 
grees of freedom (charge, spin and orbital) and the lattice 
leads to the spatial coexistence of regions with different 
properties rather than to the formation of a single ho- 
mogeneous phase. Above Tc the activated behaviour of 
the conductivitji£i£i£, the Hall conduction 8 - and the ther- 
mopowes£*&, as well as the structure of the pair distribu- 
tion function (PDF) 9 indicate the formation of small po- 
larons, i.e., of almost localised carriers within a surround- 
ing lattice distortion. Interestingly these polarons con- 
tinue to exist in the metallic phase below Tq , merely their 
volume fraction is noticeable reduced. Corresponding ev- 
idence is found in conductivity measurements^^, muon 
spin relaxationiS* 1 ^, high-resolution x-ray diffraction 14 , 
pulsed neutron diffraction 1 ^ and in x-ray absorption fine 
structure spectral. Based on these experimental obser- 
vations different scenarios for the coexistence of conduct- 
ing and insulating regions within the metallic phase of 
the manganites were discussed, which relate the metal- 
insulator transition to phase separation 4 and percolative 
phenomenaiLi&iSiSii. In particular microscopic imaging 
techniques, like scanning tunneling spectroscopji 2 ^^ or 
dark-field imaging 2 —, seem to support the latter idea. 
However, as yet the detailed nature of the spatially co- 
existing regions or phases is not known very precisely 
and even the data for the corresponding length scales is 
contradictorjiS 4 * 2 ^. 

In a recent workS& we addressed the problem of coexist- 
ing conducting and insulating regions within the metal- 
lic phase of the manganites and proposed a phenomeno- 
logical mixed-phase description, which is based on the 



competition of a polaronic insulating phase and a metal- 
lic, double-exchange driven ferromagnetic phase. Both 
phases are assumed to have an equal density of charge 
carriers and the percolative coexistence is accounted for 
by a metallic bandwidth, which depends on the volume 
fractions of the two components. The model is able to 
describe a finite polaronic volume fraction well below Tc 
and yields rather realistic x-T-phase diagrams. However, 
its sensitivity to external magnetic fields is much too 
weak and we made no attempt to describe resistivities 
or the large magnetoresistance. 



Due to the insufficient knowledge of the two differ- 
ent components of the low temperature phase, the as- 
sumptions of our previous model may be problematic, in 
particular the balance of the two components could fol- 
low from different equilibrium conditions. In the present 
work we discard the condition of equal charge density 
within the polaronic insulating and the ferromagnetic 
metallic regions and require equal pressure instead. This 
approach is well justified, since one of the components is 
insulating and the length scale of the coexistence seems 
to be short enough, as to avoid long range Coulomb ef- 
fects. We complete this new mixed-phase description by 
a model for the resistivity which is based on a random 
resistor network that accounts for the percolative nature 
of the low temperature phase. The phase diagram we ob- 
tain from the improved model is comparable to the previ- 
ous results, however, the sensitivity to external magnetic 
fields is much stronger and the ansatz for the conductiv- 
ity yields a rather large magnetoresistance close to Tc- In 
addition, the order of the phase transition from the ferro- 
magnetic metallic to the paramagnetic insulating phase 
depends on the model parameters and in particular on 
the doping x, a feature which was observed for the real 
materials 27,28 . 
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FIG. 1: Doping dependence of the polaronic binding energy: 
we compare experimental data— from conductivity measure- 
ments (dots) with the e p from our ansatz (lines). Inset: Den- 
sity of states used for the ferromagnetic metallic component. 



II. COEXISTING COMPONENTS 

In the doping range 0.2 < x < 0.5 the electronic 
properties of the manganites are dominated by the 
well known double-exchange interactiori 2 -^ 3 -^ 3 ^ and the 
electron-lattice coupling. As long as the charge carriers 
(Mn e g holes) are mobile they mediate a ferromagnetic in- 
teraction between the localised 5 = 3/2 spins formed by 
the tig electrons of the manganese. If the electron lattice 
interaction is strong enough small polarons can arise and 
the effective mass of the holes is increased by the lattice 
distortion accompanying its motion through the crystal. 
Due to this reduction of the mobility the ferromagnetic 
double-exchange may break down completely, but, as in- 
dicated by the mentioned experiments, this happens in 
a spatially inhomogeneous way. We model this feature 
by assuming a coexistence of a ferromagnetic metallic 
volume fraction, described mainly by the mean-field the- 
ory of double-exchange by Kubo and Ohata-i, and of a 
polaronic insulating volume fraction, described by an ex- 
ponentially narrow band and a paramagnetic spin back- 
ground. 



describes free fermions in an effective band of width 

W = ls [S(X + X cxt )}W , (1) 

where Wo denotes the bare band width and the field de- 
pendent prefactor is given by-i 
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To improve the above approximation of Kubo and 
Ohatai--i, we account for the orbitally anisotropic hop- 
ping 32 - 33 , which follows from the perovskite structure, 
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(4) 

and for the strong on-site Coulomb interaction. In a 
mean-field sense both is achieved 26 by working with the 
averaged density of states 



q{E)= 1 -(q + {E) + q_(E)). 



(5) 



shown in the inset of Figure ^ The densities Qq(E) be- 
long to the two bands (( — ±1), 



£k,c = -t^^cosks + Cy i^(cosfc 5 - cosfc^) 2 ) , (6) 
a 5,6' 

resulting from nearest neighbour matrix elements of 
Equation (@J with t = W/6, Equation {I}. 

We assume this non-interacting fermion model to be 
valid within the ferromagnetic volume fraction 



N 



(7) 



of the sample. Of course is temperature and dop- 
ing dependent. Its actual value is determined self- 
consistently through the equations given in Section IlIII 



B. Polaronic insulating component 



A. Ferromagnetic metallic component 



In the remaining, polaronic part of the sample, 



Within mean-field theory the quantum double- 
exchange matrix element t = t(S T + l/2)/(25 + 1) for 
the hopping of an itinerant charge carrier between two 
neighbouring sites is averaged over all amplitudes and 
directions of the total bond spin St, which is assumed to 
be placed in an inner Weiss field A = (3g[iBH* fi and an 
optional external magnetic field A cxt = f3gfj,BH^ t . Omit- 
ting orbital degrees of freedom the resulting Hamiltonian 
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= l-pO, 



(8) 



we assume all charge carriers to be self-trapped small 
polarons, i.e., their kinetic energy is exponentially sup- 
pressed and the band centre is shifted by a polaron bind- 
ing energy e p . To keep our model as simple as possible, we 
completely neglect the polaronic band width and consider 
only a dispersion- less level located at e p . In our previous 
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work 26 , using certain energy arguments, we motivated a 
doping dependence of e p of the form 



1 E x + E 2 



(9) 



Here E\ and E 2 are effective model parameters describing 
the anti-Jahn- Teller effect and the usual polaron bind- 
ing energy, respectively. If we discard the condition of 
equally charged metallic and insulating volume fractions, 
in a strict sense the old derivation is no longer valid. 
However, the comparison of the ansatz in Equation @ 
with experimental data£ for the polaronic binding energy, 
which can be extracted from the temperature and dop- 
ing dependence of the resistivity in the high-temperature 
paramagnetic phase, yields surprisingly good agreement 
over a wide doping range. Figure ^ illustrates that with 
our choice of parameters {E\, E 2 } the functional form of 
e p matches the real data quite well, as long as x is not 
too small. In absolute values the curves disagree by a 
constant factor of the order 3 to 5, which is sufficient for 
a mean-field type theory. 



III. THERMODYNAMIC STABILITY 
CONDITIONS 



where 
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denotes the spin part of the free energy (5 — S + h = 2) 
and 

vs[z] = sinh(z) coth (^) + cosh(z) (18) 

is the spin partition function. 

The set of equilibrium conditions is closed by the re- 
quirement that A is chosen such that the free energy is 
minimal. The same is required, if the above equations 
have more than a single solution. A finite value of A cor- 
responds to a ferromagnetically ordered metallic compo- 
nent, and if in addition the corresponding volume fraction 
p^' is nonzero the magnetisation of the whole sample is 
given by 



Based on the above assumptions for the two coexisting 
components of the low-temperature phase of the mangan- 
ites we are now in the position to formulate equilibrium 
conditions. The essential change to our previous work 
concerns the assumption of equal pressure for carriers in 
the metallic and the insulating region, where the pres- 
sures are obtained from 
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with (3 = l/(fcflT). For each value of the inner field A 
the equation 
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then defines the chemical potential p, which is required 
to be equal for both components. Given \i the resulting 
carrier concentrations 
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in the coexisting regions define the two volume fractions 
by the equations 

x = p U) x {f) + pWjfCiO an d p (f) + p b>) = 1 ; (15) 

which ensure the correct overall doping. Based on this 
set of equations we are able to calculate the free energy 
per site for the whole system, 

/ = X[ l - 7T cq + , (16) 



m=(l-x) p {f) SB- s [S{X + A cxt )] + P ip) SBg[SX cx 
+ x \p if) SB s [S(X + A cxt )] + p ip) SB s [SX cxt ]\ . (19) 

IV. DC CONDUCTIVITY 



The above thermodynamic relations define the phase 
boundary between the ferromagnetic metallic and the 
paramagnetic insulating phase and explain the behaviour 
of the magnetisation and of the volume fractions of the 
different components. However, they do not contain any 
information about the resistivity of the system. 

In the past, approaches which are based on the per- 
colative mixing of regions with different macroscopic re- 
sistivities have been successfully used to fit experimental 
(lat;'iil^^i. Here we follow a similar path to model the 
resistivity p of our mixed-phase system. Namely, we as- 
sume that the resistivity of the metallic component is 
proportional to the expression 
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derived by Kubo and Ohata£i, which associates p with 
the fluctuation of the double-exchange matrix element 
caused by the thermal spin disorder. The resistivity 
of the insulating component is assumed to match the 
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FIG. 2: x-T phase diagram of the mixed-phase model ob- 
tained for Wo = 2.4 eV, Ei = -0.125 eV and the two E 2 
values —0.3 eV and —0.5 eV. 



resistivity of the high-temperature phase, which in ex- 
periment is well fit by the activated hopping of small- 
polarons^i^. Hence, the resistivities of the two compo- 
nents are given by, 



P (/) = J^(ps[S(A + A ext )]+ Pmln ) 
A 
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(22) 
(23) 
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where the prefactors A and B as well as the cut-off p m i n 
arc free model parameters which could be estimated from 
experimental data. 

The resistivity of the whole sample, which may con- 
sist of an inhomogeneous mixture of both components, is 
calculated by assuming a random resistor network. More 
precisely, we choose nodes from a cubic lattice which be- 
long to the metallic component with probability p^ and 
to the polaronic component with probability p^ . Each 
of these nodes, which represent macroscopic regions of 
the sample, is connected to its neighbours with resistors 
of magnitude p«' or p^ p \ respectively. The resistivity of 
the network then yields a reasonable approximation for 
the resistivity of the inhomogeneous ferromagnetic metal- 
lic phase of the manganites. The percolative nature of 
this model, particularly in the vicinity of the phase tran- 
sition, makes the system very sensitive to small changes 
in temperature, doping and external magnetic field. 



V. RESULTS 

The numerical solution of the self-consistency equa- 
tions is rather straightforward. However, some care is 
recommended, if there are multiple solutions for the 
equations (|12|) - (|15|l . To give an example, we set the 
bare band-width and the Jahn- Teller energy equal to 



FIG. 3: Temperature dependence of the magnetisation m for 
the two parameter sets and the indicated doping levels. Nor- 
malisation: mo = S — x/2. 

W = 2.4 eV and E\ = -0.125 eV, respectively, and 
consider two typical values for the polaronic binding en- 
ergy, E2 = —0.3 eV and —0.5 eV. This choice results 
in the phase boundaries displayed in Figure [21 Without 
an external magnetic field, for each doping x, the transi- 
tion is defined by the critical temperature Tc for which 
the magnetisation m of the sample vanishes. The model 
yields reasonable values for both, the transition temper- 
atures and the critical doping x c at T = 0. 

Figure [3] shows the magnetisation m as a function of 
temperature and doping. Clearly, the order of the phase 
transition depends on both, the polaronic parameters 
{E\,E2\ and the doping x. Higher transition temper- 
atures usually correspond to a continuous, second order 
transition, whereas otherwise the transition is first or- 
der. Similar behaviour was also found for the real mate- 
rial9& 

The inhomogeneous nature of the ferromagnetic and 
the paramagnetic phase becomes evident from Figure 21 
where we show the dependence of the ferromagnetic vol- 
ume fraction pW) and of the corresponding carrier con- 
centration on temperature and doping. In the case 
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FIG. 4: Ferromagnetic volume fraction pW) (bold lines) and 
corresponding carrier concentrations (thin lines) calcu- 
lated with E 2 = -0.3 eV. 
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FIG. 5: Influence of an external magnetic field -ff^xt on the 
magnetisation m (main panel), the volume fraction p^' and 
the carrier concentration x^ (insets). 



of small Tq the sample is never completely metallic, i.e., 
a small polaronic insulating volume fraction is present 
even at lowest temperatures. On the other hand, for 
higher Tc a finite metallic volume fraction exists also 
above Tc- Since it is usually smaller than the percola- 
tion threshold the sample remains insulating. In addition 
the carrier concentration within the metallic component 
is noticeable reduced above Tc- In the case of a second 
order transition and x^ also decrease continuously 
with temperature. The presence of a metallic, double- 
exchange driven component above Tq can be related to 
the experimentally observed ferromagnetic correlations, 
sometimes interpreted as ferromagnetic clusters or mag- 
netic polarona 36 ! 37 ! 38 ! 39 . 

The sensitivity of the phase transition and of all re- 
lated quantities to an external magnetic field is illus- 
trated in Figure |3J There we show the magnetisation 
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FIG. 6: Resistivity p under the influence of an external mag- 
netic field .ffcxt for doping x = 0.3 and 0.4 and the parameter 
set {Wo, Ei,E 2 } = {2.4, -0.125, -0.3} eV. 



m together with the volume fraction p^ and the carrier 
concentration x^> of the ferromagnetic metallic compo- 
nent for the two considered parameter sets and a doping 
level of x = 0.3. Clearly, even for moderate field strength 
the critical temperature is shifted by a few degrees and 
in particular the volume fraction p^ changes noticeable 
around Tc- Of course, the latter has an important influ- 
ence on the conductivity of the system. 

Inserting the volume fractions and carrier concentra- 
tions from the mixed-phase model into the ansatz for the 
DC conductivity we obtain the resistivities p shown in 
Figure Since we are mainly interested in the general 
features of p, we set B = 1 and use the specified values 
for A and p mnl . Depending on the order of the phase 
transition p shows a sharp jump or a continuous increase 
close to Tc- This behaviour of the resistivity originates 
to a large degree from the changing volume fraction of 
the metallic component, which can cross the percolation 
threshold. However, the conductivity of the component 
itself as well as its carrier concentration strongly affect p 
for T < Tc- An external magnetic field causes a reason- 
able suppression of p, i.e., a noticeable negative magne- 
toresistance. Compared to the real compounds the cal- 
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culated effect is a bit weaker. Nevertheless, in view of the 
rather simple model for the conductivity the agreement 
is quite satisfactory. Probably, more involved assump- 
tions for the resistivities of the two different components 
could improve the magnetoresistance data. At present 
the conductivity of the metallic component is controlled 
only by spin fluctuations in the double-exchange hopping, 
which could be increased by the weak anti-ferromagnetic 
interactions active in the manganites. In addition, other 
scattering mechanisms could play a role. The adiabatic 
small-polaron approximation used for the polaronic insu- 
lating volume fraction may be questionable as well, since 
the phonon modes involved in the electron-lattice cou- 
pling do not have small enough frequency. Another po- 
tential improvement concerns the site percolation model 
used to construct the random resistor network. As was 
pointed out recently, an approach that is based on corre- 
lated percolation could be more appropriate^. An affin- 
ity to the formation of larger regions of the same type 
would naturally affect the resistivity of the system and 
its response to an external held. 



colossal magnetoresistive manganites, which is based on 
the coexistence of a double-exchange driven metallic com- 
ponent and a polaronic insulating component. Using 
modified equilibrium conditions and adding a percolativc 
ansatz for the resistivity of the mixed phase we have sub- 
stantially improved previous workSi. With realistic pa- 
rameters for the electronic band-width, the Jahn- Teller 
splitting and the polaronic binding energy our approach 
yields reasonable data for the phase boundary of the fer- 
romagnetic metallic phase and correctly predicts the ex- 
istence of a finite polaronic volume fraction well below 
the critical temperature Tc- The model shows a man- 
ifest sensitivity to external magnetic fields, including a 
large negative magnetoresistance close to Tc- 



VI. CONCLUSIONS 

In summary, we have proposed a phenomcnological 
model for the ferromagnetic metallic phase of doped 
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